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Abstract: We describe how surface interactions can affect the growth of ice crystal facets 
in contact with a substrate by lowering the normal nucleation barrier on the ice surface. We also 
describe how the resulting enhanced growth rates can produce systematic errors even when measuring 
the growth of facets not contacting the substrate. From an analysis of the diffusion dynamics we 
then develop a simple procedure for approximately correcting ice growth data, thus substantially 
reducing these systematic errors. We have found this technique to be quite useful for interpreting 
ice growth data and extracting the intrinsic attachment coefficients of ice surfaces. 

1 Introduction 

We recently developed a new experimental apparatus for making precise measurements of the growth 
rates of ice facet surfaces as a function of temperature and supersaturation [JJ [5] . During the course 
of our investigations, we found that interactions with our sapphire substrate introduced systematic 
errors that could be substantial in some circumstances. In particular, we found that ice facets that 
were in contact with the substrate tended to grow faster than otherwise expected, in a variable and 
somewhat unpredictable fashion. This enhanced growth from substrate interactions was particularly 
evident at low supersaturations. 

We believe that the enhanced growth arises because the ice/substrate interaction creates molec- 
ular steps on the intersecting ice surface, thus lowering the normal nucleation barrier that is present 
on an ice facet. A contact-angle model is sufficient to explain these observations, as illustrated in 
Figure [TJ If the ice/substrate contact angle 9 contact is less than the intersection angle 9 facet between 
the facet and the substrate, then the ice/substrate contact will act as a source of molecular steps 
on the ice facet, as seen in the left diagram in Figure [JJ The presence of these molecular steps then 
reduces the normal nucleation barrier on the ice facet, and the growth enhancement will be most 
pronounced at low supersaturation, when the normal facet growth rate is especially low. 

If ^contact > facet, however, then the contact will not be a source of steps, and there will be no 
enhanced growth arising from substrate interactions, as seen in the right diagram in Figure [TJ Since 
^contact depends on the chemical nature of the substrate, as well as on any chemical residue on an 
imperfectly cleaned substrate, we found that the ice growth behavior from substrate interactions 
can be somewhat unpredictable. 

This same effect appears to have affected other ice growth experiments reported in the literature. 
For example, in [3] the authors found that prism facets not intersecting their substrate (facet type 
P in their paper) grew approximately five times slower than facets that did intersect the substrate, 
even while the facets were all exposed to essentially the same growth conditions. Also, in [4] the 
authors found that facets in contact with the exterior of their capillary tubes grew more quickly 
than other facets, and that contact with the foreign surface also reduced the measured nucleation 
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Figure 1: An ice crystal resting on a substrate with 9 facet = 90 degrees, which is the case when 
the basal facet of a hexagonal prism rests on the substrate. At left 9 con tact < 9 facet, so substrate 
interactions introduce molecular steps on the intersecting facet surface, thus reducing the nucleation 
barrier that normally inhibits crystal growth on that surface. At right 9 contact > 9f ace t and the 
nucleation barrier is not affected. 

barrier. Qualitatively the effects described in these two experiments agreed with our observations, 
suggesting that all three experiments were witnessing essentially the same phenomenon. 

In addition, we found that substrate interactions can substantially influence the growth of neigh- 
boring ice facets that are not directly in contact with the substrate. The mechanism is that rapid 
growth of the contacting facets lowers the supersaturation field around the crystal, owing to particle 
diffusion through the background gas, and this reduces the growth rates of the neighboring facets. 
This effect is very large for background pressures near one bar, but even at 10-20 mbar (typical of 
our measurements) it can introduce significant systematic errors. In addition, as we see below, the 
effects can be important for measurements made both at low and high supersaturations. 

The remainder of this paper presents: 1) a detailed analysis of how substrate interactions affect 
the growth rates of neighboring crystal facets, those not in contact with the substrate; and 2) a 
method for modeling and approximately correcting for substrate interactions. 

2 Systematic Errors from Substrate Interactions 

Our first, rather obvious, conclusion from these observations is that growth measurements of ice 
facet surfaces in contact with a substrate may yield little useful information on the intrinsic ice 
growth rates. A faceted surface typically implies a nucleation barrier, and any reduction of that 
barrier from substrate interactions likely produces unacceptably large systematic errors. Of course, 
it is possible that substrate interactions are negligible in some circumstances, but it rests on the 
experimenter to convincingly demonstrate that this is the case. 
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We tried to find a suitable superhydrophobic coating for our sapphire substrate that would 
increase O conta ct and thus reduce the effects of substrate interactions, as shown in Figure [TJ To 
date we have not found a surface treatment that is superior to a very clean sapphire surface, but 
much research effort is currently being spent on developing superhydrophobic coatings for a range 
of materials. We believe that future measurements would likely benefit from finding suitable super- 
hydrophobic surface treatments. 

Our second conclusion is that the indirect effects on neighboring facets can be quite important 
also, even at low background pressures. Some diffusion modeling is required to understand the 
magnitude of these effects on ice crystal growth experiments. The remainder of this paper is focused 
on outlining the diffusion modeling of this problem and developing a strategy for managing and 
correcting the resulting systematic errors in experimental data. 

To develop a quantitative diffusion model, the present discussion considers the problem of crystal 
growth in the presence of an inert background gas, in which the transport of water vapor is governed 
by simple diffusion through the gas background. We further simplify to the isothermal case, which is 
a reasonably accurate approximation for growth on a substrate, since the substrate acts as a thermal 
reservoir that readily absorbs the latent heat generated from crystal growth [5]. 



2.1 The Hemispherical Case with no Substrate Interactions 

To set the notation and basic scales in the problem, we first examine the case of a hemispherical 
crystal of radius R growing on an infinite substrate. We assume a constant supersaturation 
far from the crystal and assume a single condensation coefficient a (a) for the crystal surface. The 
diffusion equation and boundary conditions for this problem are equivalent to the spherical case, 
which can be solved analytically. Following the notation in [5] , the growth rate of the crystal surface 
is 

^ — &sur fVkin® sur f 

c sa t n da , 

" ~TT7 \surf 



Csolid dR 

where a sur f = a(cr sur f) is the condensation coefficient at the ice surface and <r S urf = <j(R)- The 
diffusion equation together with these mixed boundary conditions for a (r) give the solution [5] 

_ a S urfadtff 

OLsurf + OLdiff 

where 
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where D is the diffusion constant in the background gas. For the final expression we used D a i r = 
2 x f 0" 5 m 2 / sec for the diffusion constant of water vapor in air at a pressure of one bar. To a good 
approximation D is inversely proportional to the background gas pressure. 
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Figure 2: Calculated plots of the functions a sur f(a sur f) and a meas (<7oo) denned in the text, for 
the case of a hemispherical crystal growing on a substrate. The individual points show the corre- 
spondence between the two curves. Here we used otdiff — 0.25, corresponding to a crystal radius 
of R as 30 /zm in a background pressure of 20 mbar (both typical of our measurements), and an 
intrinsic attachment coefficient of a sur f (a sur f) — exp(— 0.015/<r sur /) (typical for basal facet growth 
at temperatures near -13 C) [2]. 



From the analytic solution we also have the supersaturation at the ice surface 

Otdiff 



&surf 



Otsurf + Otdiff 



and eliminating a sur f from this expression gives 



@ surf ^"o 



OtdiffVkin 



(2) 



(3) 



which is exact for the spherical and hemispherical cases. 

When measuring ice crystal growth rates [1] [2] , we typically establish a supersaturation cr^ far 
from the crystal and then measure the crystal dimensions as a function of time and <7oo , with the 
goal of extracting the intrinsic attachment coefficient a {(J surf) for the growing surfaces. In the limit 
of fast diffusion (or slow attachment kinetics) , a sur f <C otdi / / and the hemispherical case gives 



Otsurf ~ Oi meas 



V kin® og 



(4) 
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which defines the measured quantity a meas . If diffusion cannot be neglected, then 

_ Ot sur f Ct<iiff 

Corneas — ' \&) 

OLsurf + Otdiff 

and a plot of a meas (ctqo) typically deviates from the desired a sur f (a sur f) at high (Too. An example 
of this is shown in Figure [2J and plots of a meas (jJoo) from our data show this same functional form 

mm- 

The example in Figure [2] shows that a m eas (coo) is approximately equal to a SU rf (cr sur f) at 
low supersaturations, which is precisely when a SU rf <C adiff- Increasing Coo increases a SU rf until 
Q-surf (&surf) becomes larger and the criterion a sur f <C otdiff is no longer satisfied. When this 
occurs, the two curves in Figure [5] deviate from one another. 



2.1.1 Correcting the Hemispherical Case 

In our experiments we determine and measure the growth velocity v of a facet not contacting 
the substrate, and from Equation @] we derive a meas (o"oo) entirely from measured quantities. Our 
goal, however, is to determine the intrinsic condensation coefficient a sur f (a sur f) . For our idealized 
hemispherical crystal, we could apply a correction to the data that transforms a meas (coo) into 
((Jsurf) using 

1 Coo (b) 

a di ff J 

/ \ / 1 Corneas \ 

Crsurf \Gsurf ) - — I -L I Q-meas 

V a diff J 

where adiff is derived from the measured crystal radius using Equation [1] Clearly this correction 
works best when a meas "C adiff, or equivalently a sur f <§; adiff- As long as the crystal growth 
is limited mainly by attachment kinetics and not by diffusion, then we can use this correction to 
reasonably extract a sur f (a sur f) from the measured a meas (coo). 



2.2 Plate Growth with Substrate Interactions 

We next consider the growth of simple hexagonal plate crystals, which is what we observed in [2], In 
this case one basal facet is resting on the substrate, and we measure the growth of the opposite basal 
facet along with the growth of the six prism facets. We typically have Vbasai ^ v pr i sm , and there is 
no analytic solution of the diffusion equation that gives the supersaturation field cr{x) around the 
crystal. 

To see the problems that result, consider the example shown in Figure |31 in which we examine 
the basal growth of a plate of radius R 30 /im and thickness H « 4 /im, again at a pressure 
of 20 mbar. If it so happens that a pr i sm (a) = a^asai (&) , then v pr i sm ~ Vbasai and we see that 
the resulting a rneas (a^, ) looks quite similar to the hemispherical case in Figure [21 with some adi / / 
appropriate for the crystal size and geometry. In particular, we see that a meas (coo) ~ a sur f ((J SU rf) 
at low supersaturations, as we saw with the hemispherical case. 

If we increase the prism growth rate, however, in this example by setting a pr i sm — 1, then the 
faster prism growth pulls down the supersaturation around the entire crystal. Because a^asai (c) is 
a strong function of er, this means that the measured basal growth rates are substantially reduced. 
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Figure 3: Calculated plots of the functions a sur f(a sur f) and a m eas (coo) for growth of the top basal 
surface of a hexagonal plate crystal, with the lower basal surface resting on the substrate. The indi- 
vidual points show the correspondence between the three curves. In this example we used the model 
of a circular disk crystal of radius R « 30 fim and thickness H w 4 /jm in a background pressure of 20 
mbar, with the intrinsic basal attachment coefficient a basa i = a sur f (cr sur f) = exp(— 0.015/cr stir y). 
For the middle curve we used a prism = a basaU which gives v prism w v basa i and an a meas (cx^) 
curve that is similar to that shown in Figure [5J For the lower curve we used a pr i srn = 1 , giving 

Vprism ^ ^basal and an CXrneas 

((Too) curve that is substantially distorted, especially at low a^. 

Furthermore, the reduction factor is especially large at low supersaturations, as can be seen by 
examining the curves in Figure [3J Now this is a rather extreme case, with exceptionally fast prism 
growth, but it demonstrates how the measured a rneas (coo) can be distorted at both high and low 
by diffusion effects, even at quite low background air pressures. These effects can make it difficult 
to reliably determine a sur f (cr sur f) from the measurements, requiring diffusion modeling. 

2.2.1 Correcting Experimental Data 

In principle one could use numerical modeling of particle diffusion to extract a sur f (cr sur f) from 
ct-meas {? oo) , provided the difference is not too large. In practice, however, such modeling is quite 
computationally intensive, and it is problematic to do a detailed numerical model for each of hundreds 
of test crystals. We would therefore like a simpler analytic approach that would allow us to apply a 
rough correction to our data. To produce such an approximate correction, we generalize Equation 
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[3] by writing 

~ \ ' fiVi (7 s 

&sur f ~ &corr — ^oo / \ > ) 

~^ adiff,effVkin 

where the sum is over the different growing regions of the crystal (in this case the one basal facet 
and six prism facets). The fi give the fractional areas of the different regions (with ^2 fi = 1)> where 
each region is growing with velocity Vi. From Equation [T] we take 

0.15 (^) (£-) 

where 

i 

and the Ai give the areas of each crystal region. Note that Equation [7] is exact for the hemispherical 
case, where it reduces to Equation [3] For the facet of interest we then have 

O^surf \&surf ) ^ &meas 
VkinCsurf O corr 

and this correction can be applied to each data point. 

A key feature in this correction is that the Vi, A i: and /, can all be extracted from the measured 
crystal size and geometry as a function of time. If the prism growth is influenced by substrate 
interactions, as we describe above, we can still use this formalism to correct for changes in the 
supersaturation field that affect the basal growth measurements. As long as the difference — a corr 
is not too great, we can determine a sur f (a sur f) for the basal surface. 

Figure [4] shows how this correction process works in practice, using a plate-like crystal growing 
at a temperature of -12 C, taken from [2]. Note that the prism growth rate at the beginning of 
this measurement was v pr ism ~ 63 nm/sec, when doo sa 0.33%. If a pr i S m had been equal to unity, 
as we assumed in the example in Figure [31 then the radial growth would have been a much faster 
Vprism ~ 900 nm/sec. Thus we verify that the a pr i sm — 1 case in Figure [3] was somewhat extreme, 
and in our actual data the corrections are usually substantially smaller. 

Note also that our fit to the uncorrected data in Figure @] gave Q-diff = 0.09, a value roughly 
a factor of two smaller than the adiff calculated from the size and geometry of the crystal. We 
observed this discrepancy frequently in our basal growth data, and we find it is nicely resolved by 
applying the above data correction. The corrected data give a significantly better measure of the 
intrinsic attachment coefficient a sur f (<J SU rf), with reduced systematic effects compared with using 
only the uncorrected a me as (coo). This improvement is accompanied by only a modest increase in 
the scatter of the data points arising from the imperfect correction process. 

It is straightforward as well to apply this same correction process to measurements of prism facet 
growth, where one prism facet rests on the substrate. The area factors fi must be changed to fit the 
prism growth case, but otherwise the prism growth correction is identical to the basal case described 
above. 

3 Conclusion 

The approximate correction described above serves several purposes in the analysis of our data: 1) 
it reduces the systematic errors in our determination of a SU rf {(J surf) for the growth of facets that 
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do not contact the substrate; 2) it eliminates the fit parameter otdiff when fitting the observed 
o^meas (c<x>) (see [2]); 3) a comparison of a meas ((Too) and a SU rf (a surf) indicates which crystals have 
large correction factors and thus should be given less weight in determining a sur f (a sur f) ; and 4) 
the correction allows us to make reliable measurements of a sur f (a sur f) over a broader range of 
experimental conditions than we could without the data correction. 

Although the correction is only approximately valid, it nevertheless does an adequate job of 
reducing an important systematic effect. We have found it to be a valuable tool for interpreting 
measurements of the growth of ice crystals on substrates at low background pressures. 
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Figure 4: An example of actual ice crystal growth data corrected using the procedure de- 
scribed in the text. (Top) The measured radius and thickness of a hexagonal plate crystal as 
a function of time during a growth run from [2], along with 1500(700 (in absolute units). The 
background pressure during this measurement was 25 mbar. Note that at early times the ra- 
dius increased more rapidly than the thickness (v pr i sm > Vbasai) because of substrate interac- 
tions. (Middle, Bottom) The directly measured condensation coefficient a meas {a oo) (open sym- 
bols) and the corrected measurement a(a sur f )g (closed symbols). The fit lines in both these 
plots are a rneas (a 00 ) = exp(-0.0205/cr)a di ///(exp(-0.0205/cr) + ay*//), with otdiff = 0.09, and 
Oi{cr stlr f) = exp(— 0.0185/er). 



